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We introduce a new representation of an integer spin S via bosonic operators which is useful 
in describing the paramagnetic phase and transitions to magnetically ordered phases in magnetic 
systems with large single-ion easy-plane anisotropy D. Considering the exchange interaction be- 
tween spins as a perturbation and using the diagram technique we derive the elementary excitation 
spectrum and the ground state energy in the third order of the perturbation theory. In the special 
case of S = 1 we obtain these expressions also using simpler spin representations some of which 
were introduced before. Comparison with results of previous numerical studies of 2D systems with 
5* = 1 demonstrates that our approach works better than other analytical methods applied before 
for such systems. We apply our results for the elementary excitation spectrum analysis obtained 
experimentally in NiCl2-4SC(NH2)2 (DTN). It is demonstrated that a set of model parameters (ex- 
change constants and D) which has been used for DTN so far describes badly the experimentally 
obtained spectrum. A new set of parameters is proposed using which we fit the spectrum and values 
of two critical fields of DTN. 
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I. INTRODUCTION 



Many analytical approaches to consideration of quantum magnetic systems with localized spins base on repre- 
sentations of spins via bosonic or fermionic operators allowing to turn from spin Hamiltonians to those describing 
ensembles of bosons or fermions.^ At low temperature the form of the representation depends on the ground state of 
the spin system under discussion. The Holstein-Primakoff and the Dyson-Maleev representations are frequently used 
for magnetically ordered phases, the Jordan- Wigner transformation proved to be useful for S — 1/2 spin chains, the 
bond-operator formalisnP was proposed for spin liquids with singlet ground state, etc. 

We propose in the present paper a new representation of an integer spin S. This representation should be useful 
in describing the paramagnetic phase in which all spins are mainly in the quantum state with zero quantum number 
for projection onto a preferential direction. As an example of particular system of this type we choose a Heisenberg 
magnet with large single-ion easy-plane anisotropy which Hamiltonian has the form 

^ = ^E-^^^S^S, +i?^(5ff , (1) 

where summations are taken over all sites of the lattice with arbitrary spatial dimension, D > and signs of exchange 
constants J^- are not important in the following. The first term in ([ij is supposed to be small enough compared 
to the second one so that the system to be in the paramagnetic phase at T = 0. There is a quantum critical point 
(QCP) D = Dc separating the paramagnetic phase at D > Dc and a magnetically ordered or a spin-liquid one at 
D < Dc- The value of Dc as well as the phase type at D < Dc depend on the lattice spatial dimension and the 
exchange interactiorPI^. 

Quite a few compounds can be mentioned with the paramagnetic ground state of the considered type which are 
described by the Hamihonian ^: CsFeBr^tS!^ CsFeClJ^, SraNiPtOj^, NENCpHISl^ NENPESEl, NBYCP, NiCla- 



4SC (NH2 ) jj^itSSl (a pparen tly the most intensively studied compound of this type in recent years), (Ni(C5H5NO)6)(N03)2 
and NiSnClg • 6H2CP2I2II. AH the mentioned materials have 5* = 1 and all of them are quasi-lD magnets. It is prob- 
ably the reason why the majority of the recent theoretical investigations of the model ([T]) with large D > focus 
on weakly coupled or independent spin chains'^ In particular, the elementary excitation spectrum has been 

derived before using a random phase approximatiorPl^, the regular perturbation theorj^^H^ a "generalized spin- wave 
approach"-^ and some other self-consistent procedureJ^ 

Treating the first term in Eq. ([T]) as a perturbation, using the proposed spin representation and the conventional 
diagram technique we derive below expressions for the elementary excitation spectrum and the ground state energy 
in the paramagnetic phase in the third order in the perturbation theory. For brevity, this approach is referred to 
hereafter as an expansion in terms of J/ D while one should bare in mind that constants Jtj in Eq. ([I]) are assumed 
to be arbitrary. We also obtain these results in the special case oi S = 1 using simpler spin representations some 
of which were introduced before in Refs ;'^^'^^ . In the particular case of a spin chain with S = 1 our expression for 
the spectrum coincides with that obtained in Ref.f^Sl only fo r this special case using the regular perturbation theory. 
Comparing our results with those of numerical calculation^^^ which were carried out for S = 1, square lattice 
and antiferrom agnetic exchange we show that our approach works better than other theoretica l met hods proposed 
gQ faj ^5 | 38 | 4 0i44j ^ particular, our results are in very good agreement with the numerical one^^^^^ not very close 
to the QCP. At D Dc only the spectrum of long-wavelength excitations is reproduced unsatisfactorily that is a 
consequence of strong fluctuations near the QCP. 

We demonstrate that our approach is applicable to the intensively studied compound NiCl2-4SC(NH2)2 (DTN^^^^H 
described by the model (fij). It is shown that our expression for the spectrum describes badly the spectrum obtained 
in the neutron experimenP^ with the conventional set of parameters (values of exchange constants and D) attributed 
to DTN before. A new set of parameters is prop osed for DTN using which we fit well the neutron data and reproduce 
values of critical fields found experimentally221, should be noted that in comparison with the conventional model 
used for DTN discussion we take into account also an exchange interaction between two magnetic sublattices of DTN 4^ 
which should be taken into consideration as the recent ESR experiment^® demonstrates. 

The rest of the present paper is organized as follows. In Sec. |ll] the new representation of an integer spin is 
proposed and details of the diagram technique based on this representation are discussed. In Sec. |III| expressions for 



the elementary excitation spectrum and the ground state energy are derived in the third order in J/D. In Sec. IV 
we compare our results with those obtained before using other approaches. We discuss also representations which 
are val id o nly for S = 1 some of which are simpler than the general one. The expression for the spectrum derived 



in Sec. Ill is applied for analysis of the DTN spectrum obtained in the neutron experiment!^. Sec. |v] contains our 



conclusion. 
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II. REPRESENTATION OF AN INTEGER SPIN AND TECHNIQUE 

The ground state of a system described by the Hamihonian ([T]) in the hmit of J/Z3 —> is a direct product of 
states \Sf = 0): 11^ (g) | S'f = 0) . The lowest excited states can be constructed from the ground state by substituting 
\Sf = ±1) for \Sf = 0) at any i. The energy of such states is equal to D and one leads to the doubly degenerate 
dispersionless elementary excitation spectrum 

eop = D. (2) 

The exchange interaction gives rise to the spectrum dispersion. When it is small enough one can find expressions 
for the spectrum and other observables considering the exchange as a perturbation. In particular, the elementary 
excitation spectrum of the spin chain with 5* = 1 and the exchange coupling between only nearest neighbors was 
calculated in Ref.l^ in the third order of the regular (non diagrammatic) perturbation theory. In contrast, our aim 
is to construct a spin representation which opens the door to calculations utilizing the diagram technique and which 
allows to obtain such expressions for observables easier. In particular, we recover below the result of Ref.'^^ using the 
diagram technique. 

We propose the following expressions for projections of an integer spin S: 

S- = nts - na,i, (3) 



y 1 + "b,i y 1 + ^^a,^ 

where and bi are bosonic operators, n^^i = a\ai^ rih i — blbi and 

{S^ = -{n + 1)1 4 \St = -n) = Vr^, 
{S^ = n + 1| 6j ^ n) ^ V^TTT, 
a, \Sl=Q) =0, 

n > 0. 

Operators a) and W create excitations with = —1 and +1, respectively. The subspace of physical states is 
constrained by the following conditions: 

na,inb,i = 0, (5) 
nt. < S. 

It can be readily verified that representation Q-Q reproduces the spin commutation relations on the physical 
subspace defined by Eqs. ^ and 

Condition ([s]) selects states having at any site only excitations of a- or b- type. This constraint can be satisfied by 
adding to the Hamiltonian the term describing an infinite repulsion between a and b particles 

= ^ ^ 4bla^b^, U ^ +w. (7) 

i 

After this modification matrix elements of operators ([Sj-Q become zero between states from physical and unphysical 
subspaces constrained by conditions ([6|. It means that at zero (and most probably at low) temperature we can use 
Eqs. (|3])-(4| and the diagram technique forgetting about condition (|6|, as it is done in a similar situation in the case 
of the Ho stcin-Primakoff representation. We prove this statement below for = 1 by performing calculations using 
Eqs. ([3|-(4]), taking into account constraint term Q and introducing to the Hamiltonian the additional term 

= ■^^(aI4a.a. + 5l6j6,6,) , C/ ^ +cx), (8) 

i 

which explicitly selects states with no more than one a or b particles as condition ^ requires at = 1. It is worth 
noting that we could construct a spin representation similar to (|3|-(|4]), which matrix elements are zero on states from 
physical and unphysical subspaces so that it was not necessary to introduce term Q to the Hamiltonian. However 



4 



such a representation would be very cumbersome. On the other hand term ([t]) does not complicate calculations much 
at any integer S. That is why we use below Eqs. (|3])-(|4| with constraint term ([t]). 

At sufficiently small exchange interaction and low temperature we expect densities of a and b particles to be small. 
Therefore one can expand square roots in Eq. Q into series and restrict oneself by the first terms of the resultant 
normally ordered expressions which have the form 



5+ bl (^ci - C2 bjbi^ + (ci - C2 aja^^ a^, 
ci = ^S{S+1), 

Using Eqs. ([s]) and ([9| and taking into account Eq. ^ we obtain from Eq. ([T]) 
-H = ^ eip {alap + 6^6p) 



+Ei"^pK^-p+«p^-p 



(9) 
(10) 

(11) 

(12a) 
(12b) 



p 
1 

N 
1 

N 



E 

?2=PS 

E 



Pl+P2=P3+P4 



D 



ClC2 



( Jl + ^3) 



(^a\ala3a4^ + &| 62^3^4) + [U - J3-1] al4«3&4 I" (12c) 



C1C2 



Ji 



where Jr, 



Pl+P2+P3=P4 

J2j JijC^^^^^ , N is the number of unit cells, 



(^b\a\ala4 + a\b\blb4 + a\a3a2bi + b\b3b2ai^ , 



(12d) 



eip = eop + ^Jp = D+^Jp, (13) 

eop is defined by Eq. ^ and we omit indexes p in Eqs. (12c I and (12d). Note that we take into account only terms 
with no more than four operators in Eq. ( |12[ ) . It can be shown that terms containing more than four operators which 
appear from higher order terms in the series expansion of square roots in Eq. Q lead to contributions to the spectrum 
and to the ground state energy of the order of {J/D)^ and higher. As our aim is to calculate these quantities only up 
to the third order in J/D, we can use Hamiltonian (12). 

It is convenient to introduce the following Green's functions: 



Gip) = ~i{apal) = -i{bpbl), 
Fip) = -i{blj,al), 



(14) 
(15) 



where p = (w,p) and Up is the Fourier transform of ap{t). Naturally, the equality (opaj,) = {bpblp) is satisfied. Dyson 
equations for these Green's functions have the form 



G{p) = Go(p) [1 + EpG(p) + UpFip)] , 
Fip) = G„{-p) [npG(p) + S-pF(p)] , 



(16) 



where Go{p) = {lo — eip + i6) ^, Ep and Hp are normal and anomalous self-energy parts, respectively. The solution 
of Eq. ( 16 1 has the form 



G{p) 

F{p) 
V{p) 



-ip 



_V{p) ' 

— Cip — eip(Ep + S_p) + a;(I]_p — Ep) — EpE_p 



in,„ 



(17) 

(18) 
(19) 



III. APPLICATION OF THE APPROACH 



In this section we apply the method described above for calculation of the elementary excitation spectrum and the 
ground state energy. 
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FIG. 1. (a) Diagrams for the normal self-energy part Ep giving non-zero contributions of the second order in J/D. Solid and 
dashed lines stand for Green's functions G{p) of a and b particles, respectively, defined by Eq. (171. Lines containing solid and 
dashed parts stand for anomalous Green's functions F{p) defined by Eq. ( |18[ ). Bare vertexes are defined by Eqs. ( |12c[ ) and 
(12dl. Open and filled circles represent renormalized vertexes equations for which are presented in panel (b). 



A. Elementary excitation spectrum 



The elementary excitation spectrum ep is defined by poles of Green's functions (17) and (18) 

I?(ep,p) = 0. 



(20) 



Let us consider diagrams for the normal self-energy part Ep some of which are shown in Fig. ^ If a diagram contains 
a contour that can be walked around while moving by arrows of functions Go{p), integrals over frequencies in such a 
diagram give zero!^ In diagrams without such contours there are at least two vertexes ( |12d[ ) or at least one vertex 
(12d) and one anomalous Green's functions F{p). As the vertex (12d) is of the order of J/D and the numerator of 
F{p) is 0{J/D) (see Eq. (12b)), the normal self-energy part is the value of the order of [J/DY. Thus, we obtain for 
Ep in the first order in J pJ 



0. 



(21) 



Let us turn to diagrams for the anomalous self-energy parts Hp. One of them is shown in Fig.[2|^b). The contribution 
to Hp of the first order in J/D is given by term ( 12b). Diagrams for Hp with one vertex contain sums like Jp+k — 
in the first order in J/D. Therefore such diagrams are at least of the second order in J/D. A consideration of other 
diagrams similar to that presented above for Ep shows that they give contributions of the second order in J/D and 
higher. Then we have in the first order in J/D 



It is convenient to rewrite Eq. (20 1 using Eq. (19) in the form 

9 / , ^ n2 1^ |2 



It follows from the previous discussion that 

Cp - eip 



ep - eip - E(,p^p)) (E(,p^p) - E(_,p^_p)) 



5^(.p,P) = 0{{J/Df), 



(22) 

(23) 
(24) 



S(ep,p) - S(_,^,_p) = 0{{J/Df). 



Bearing in mind these equations and that our aim is to find the spectrum in the third order in J/D we can write 
Eq. (l23l) in the form 



n, 



(ep,P) I 



0{{J/Df). 



(25) 
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As follows from this equation, the first corrections to eip are of the second order in J / D. Taking into account that the 
first corrections to self-energy parts are also of the second order in J/-D, we can replace (ep,p) by (eip,p) in Eq. (25) 
and write down the final formula for the spectrum calculation up to the third order in J / D 



Let us find the spectrum in the second order in J/D. As is clear from Eq. (261, one has to use the anomalous 



(26) 



self-energy part in the first order in J/ D for which we obtain Eq. \22\ . The normal self-energy part has to be found in 
the second order in J/D in which diagrams presented in Fig.jlja) should be taken into account. In order to calculate 
these diagrams it is necessary to find vertexes denoted by open and filled circles in the zeroth order. As follows from 
the diagram analysis presented above, the zeroth and the first order contributions to the vertexes are presented by a 
series of ladder diagrams only (see Fig. [ijb)). Diagrammatic equations in Fig. [l|b) have the following explicit form 
in the zeroth order in J/D: 



r(") = D-2 



2D-n' 



r(o) 

^ 2 

2D- 



u 



(27a) 



(27b) 



where f2 is the sum of incoming lines frequencies. The solution of Eq. (27) at C7 — >■ +oo has a simple form 

D{2D 



.(0) 



n) 



.(0) 



AD-n 
2D-n. 



(28a) 
(28b) 

After substitution of Eqs. (28 1 to diagrams for self-energy parts and integration over frequencies f2 gives a difference 
of the kind eik^ — cika with some momenta ki^2- Then corresponds to values of the order of J/D in diagrams for 
self-energy parts and it can be neglected in the vertex calculation in the zeroth order in J/D. As a result we obtain 
for vertexes from Eqs. (28) 



ri°^ ^ D/2, 



AO) 



2D. 



(29a) 
(29b) 

One obtains for the normal self-energy part in the second order in J/D as a result of calculation of diagrams shown 
in Fig. [ija) taking into account Eqs. (29) 

2cf + iCiC2 - 1 ^ ^2 



;;(2) 

^(eip,P) 



8D 
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Using Eqs. (30), (26 1 and (22) we have for the spectrum in the second order in J/D (cf. Eq 

2c\ 



e2p 



.(2) 



1 



Op 



n 



(1) 

(eip.p) 



^ D 



4ciC2 

81) 



_2 _ \ ' 72 _ _ 



P 



(31) 



In order to calculate the spectrum in the third order in J/D one has to find the normal and anomalous self-energy 

parts in the third and the second order, respectively. To find Ep in the third order we have to take into account 

diagrams shown in Fig. l2ja). Besides, one has to consider also diagrams presented in FigJIJa) taking into account 
. „ . . ........ . . 



the first order corrections to the vertexes for which we have after solving equations in Fig. fTTb) 



(1) 



^'-D-'-n 

2 8 



1 



:J: 



C1C2 



3-1 



( Jl + J3), 



(32) 
(33) 

The only diagram of the second order in J/D for the anomalous self-energy part is shown in Fig. [2]jb). It should be 
found using Eq. ( 33 1 . As a result of straightforward calculation we obtain 



= 2D-n- J- 



3-1- 



.(3) 

^(eip.P) 



16D^ 



J^2 ^ JqJkJq-k 
k,q 

I ~ 2clc2 



4ciC2 - C^C2 

8D^ 



8D^ 



2l^J-Sr T T T 

^2 -'p-'k-'q- 

k,q 



k+p, 



n 



(2) 

(eip.P) 



AD N 



^5I^kJk-p + |^^$^Jk- 



(34) 



(35) 



Using Eqs. (26), (21|, (22|, (30), (35) and (34), one has for the spectrum in the third order in J/D 

1 



e3p 



D 



+ 



Ci(ci +C2) 

i6i:i2 



2cf + 4ciC2 - cfc^ 1 



8D 



72 _ 



N ^ 

k 



r 

SDN P 



1 



+ 



3c^C2 



^ ^ t/q(/kt/q_k 
k,q 

2CiC2+ClC2 1 



4ct 



4ciC2 
81)2 



8D2 



N2 ^ 
k,q 



cf 1 

>/p>/k>/q-k+p — g£)2 '^P ]y J]iJ\i-p 



t3 

16L>2 -"p- 



(36) 



The spectrum has a minimum in some point p = Po, in which it is separated from the ground state by a gap. 
For example, this minimum is located at p equal to the antiferromagnetic vector po in the simple square or the simple 
cubic lattices with antiferromagnetic exchange interaction between only nearest neighbors (po = (7i',7r) and (7r,7r,7r), 
respectively). The gap decreases with decreasing D and it vanishes at QCP D = Dc- Upon further decreasing of 
D a "condensation" takes place of elementary excitations with momentum pg which corresponds to appearance of 
long-range magnetic order. We show below by particular examples that aX D > Dc series in J/D for some quantities 
converge slowly. It is a manifestation of strong fluctuations near QCP. 

It is worth noting that Ci, C2 ~ 5 at 5 1 (see Eqs. (10) and (111). Therefore, as follows from Eq. (36 1, the 
expansion parameter is actually S^J/D (not J/D) at S* ^ 1. It means, in particular, that D^ ^ S^J in the case 
of antiferromagnetic exchange interaction between nearest neighbor spins on the simple square or the simple cubic 
lattices. This conclusion is consistent with the result of the spin wave analysis carried out in antiferromagnetic phase 
(see, e.g., Ref.l^. Thus, a very strong anisotropy or a very small exchange are required for the paramagnetic phase 
stability at 5* ^ 1 and the paramagnetic phase is absent in the classical spin limit at any finite J and D as it must 
be. 



B. Ground state energy 



It is useful to calculate the ground state energy for the sake of comparison with numerical results. Diagrams giving 
non-zero contributions of the second and the third order in J/ D are shown in Fig. [3j The straightforward calculation 
of these diagrams leads to the following expression for the ground state energy: 

^''^"SDN^'^P' lefe ^ ^P^k Jp+k. (37) 

P P,k 





FIG. 3. Diagrams giving non-zero contributions of the second and the third order in J /D to the ground state energy. Same 
notation as in Fig. [l] 

Notice that the first non-zero term in Egg is of the second order in J/D. 

IV. DISCUSSION AND COMPARISON WITH PREVIOUS RESULTS AND EXPERIMENT 

Note that we do not specify the type of exchange interaction Jij in Eq. ([l]) and the lattice type and dimension 



while deriving Eq. (36). Then Eq. (36) is applicable, in particular, for the spin chain with S = 1 and with exchange 
interaction between nearest neighbors only. The spectrum in this special case was calculated before in the third order 
in J/D in Ref."^^ using the regular perturbation theory. It is easy to verify that the result of Ref.^^ coincides with 
Eq. (p6|) in this case. 



A. Other spin representations for S = 1 



In the particular case of 5 = 1 a number of simpler spin representations can be introduced. In one of them Sf is 
^fiven by Eq. ^ and Sf has the form (cf. Eq. (|4|) 



S+ ^V2 b. 



(38) 



In contrast to Eqs. ([3])-(|4| the introduction of term ([s]) into the Hamiltonian is necessary now because Eq. (38) has 
non-zero matrix elements between states from physical and unphysical subspaces. Representation (38) is actually 
introduced in Ref.'^, where spin chains in magnetic field are discussed. We have made calculations of the spectrum 
and the ground state energy using this representation and recovered Eqs. (36) and (37) at S = 1. 

It should be noted that at 5 = 1 one can make calculations without taking into account the constraint term ([s]) by 
adding to Eq. (38) projector operators 1 — na,i — Ub^i as follows: 



V2 (1 



(39) 



In contrast to Eq. ( 38 ) this representation has zero matrix elements between states from physical and unphysical 
subspaces and there is no need to take into account Eq. (|8| now. We have carried out calculations with Eq. ( 39 ) and 
recovered Eqs. (p6|) and ([37| at S = 1 once again. 

1 is introduced in Ref.l^ 



The following spin representation for S ■ 



nbA h + ajyr 



(40) 



which is equivalent to our one (39) on the physical subspace and which, therefore, should lead to the same results. But 
authors of Ref.l^ do not take advantage of the opportunity to find physical quantities in the form of series in terms of 
powers of J/D which such representations provide. They expand roots in Eq. (40) and analyze the spectrum of the 
Hamiltonian in the harmonic approximation. They also find renormalization of this spectrum by taking into account 
the simplest diagrams and carrying out some self-consistent calculations for the square lattice and antiferromagnetic 
exchange. We compare below results obtained in Ref.l^ with our ones. 



B. Other approaches to the spectrum calculation 



After substitution of the self-energy parts (21 ) an d \22\ o btained in the first order in J/D to the general expression 
(23) for the spectrum we recover Lindgard's formulaf^'^^^ 



eip = + 5(5+1) Jp], 



(41) 
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which was found in Refs!^'^ within the random phase approximation. But as it is clear from the above discussion, 



Eq. (411 is correct only in the first order in J/D. The second and the third order terms in J/D in Eq. (41) differ 



significantly from those in Eq. (36). Case study that is done below shows that Eq. (41) works quite badly when J/D 
is not very small. 

A "generalized spin- wave approach" (GSWA) is used in Refs.E^f^H for the spectrum consideration in DTN. In the 
framework of this approach in which fluctuations are taken into account in a mean field fashion the following expression 
is obtained: 



(42) 



where parameters s and are determined as a result of self-consistent calculations using equations 



D 



-y 



Jo 



-y 



p 



gswa 



gswa 



(43) 
(44) 



It is seen that E q. (|42| is a modification of Eq. ( [4T| ) at 5 = 1: D and Jp are renormalized. We show below by case 
study that Eqs. (42 )-(^44p work much better than Eq. (41) but Eq. (36) proves to be more precise and convenient 



because the numerical solving of integral equations and numerical calculations are not required. 

It should be noted that an a pproa ch very similar to that leading to Eqs. (42 )-(44) is proposed in Refs.'^^'^. However 
the resultant equations in Refs.'^^'^ are more cumbersome than Eqs. (42 )-(|44) whereas they work slightly worse than 
Eqs. (42 )-([44|) as our comparison shows with the numerical results of Refs.lSsEZl, That is why we do not consider here 
in detail results of Refs.'^^ 



Comparison with numerical results 



The elementary excitation spectrum is found in Ref.^^by Monte-Carlo calculations for the square lattice and S* = 1. 
The exchange i ntera ction is taken in Ref.l^ to be positive (antiferromagnetic) and equal for all nearest neighbors. 



1 in the whole Brillouin zone when D is not 
The phase with the long-range magnetic 
1361) describes well the 



Results of Refs.'^^'^ together with the spectrum calculated using Eqs. (36), ([42|-(44) and (41) are shown in Fig. |4] 
for D = 6J and D — 10 J. It is seen from Fig. [4] that Eq. (36) works we ' ' ------- - - - 

very close to the critical value D^^ « 5.65 found numerically in Refs.'^'^ 
order (antiferromagnetic phase) is stable at D < D^^ . Notice also that even at D > Eq. 
spectrum of short-wavelength elementary excitations. 

The noticeable deviation of the long-wavelength excitation spectrum given by Eq. p6[ ) (with momenta p w po = 
(tTjTt)) at > is a result of strong fluctuations near QCP which manifest themselves in a bad convergence of 
the series in terms of powers of J/D. It is illustrated by Fig. [5] which shows the dependence of the spectrum gap at 
p = Po on D. It is seen that the second and the third order corrections to the gap are approximately equal to each 
other at Z? > D^^ and the v alue of D^ found from the relation e^p^ = is 12% smaller than the value of D^^ w 5.65 
obtained numerically in Refs.'^'^^^^. 

The result of the ground state energy calculation is shown in Fig. |6j which demonstrates that Eq. (37) works well 
when D is not very close to D^^ and that series for Egs converges slowly at Z? > D^^ . 

Comparison of the results obtained within different approaches and presented in Figs. |4]-[6] shows that formulas 
derived above using the method suggested in the present paper, on the whole, work better. At the same time it 
should be noted that despite its simplicity GSWA gives quite precise results. 



D. Application to NiCl2-4SC(NH2)2 

At the present time the most extensively studied compound described by the model ([T]) and having the paramagnetic 
ground state is NiCl2-4SC(NH2)2 which is known as DTlNpSHIII, xhe magnetic subsystem of DTN consists of Ni ions 
with S = 1 and the Lande factor g — 2.26. Magnetic ions form a body-centered tetragonal lattice which can be 
viewed as two interpenetrating tetragonal sublattices. The exchange interaction between spins inside one sublattice 
is antiferromagnetic and strongly anisotropic: the exchange constant along the tetragonal hard axis (z axis) is much 
larger than those along x and y axes. That is why DTN is considered as a quasi-lD compound. Hamiltonian M with 
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FIG. 4. (Color online.) Elementary excitation spectrum ep of the model ([T| on the square lattice with S = 1 and antiferromag- 
netic exchange J > between nearest neighbors for (a) D — & J and (b) D — lOJ. Monte-Carlo data of Refs.'^^'^ are shown 
by points. Dash, dash-dotted and solid lines are drawn using Lindgard's formula (41 1, GSWA formula (42 1 and Eq. (36 1 (the 
present paper result), respectively. The deviation of the spectrum ( |36[ ) from the numerical data at D = 6J in the close vicinity 
of the minimum at (tt, tt) is a consequence of the proximity to the QCP D = D^^ ~ 5.65 separating the paramagnetic and the 
antiferromagnetic phases (see the main text). 




FIG. 5. (Color online.) The spectrum gap at p = (tt, tt) is presented as a function of D in the model ([TJ with S = 1 and 
antiferromagnetic exchange interaction J between only nearest neighbors. Monte-Carlo data of 

Ref.™ 

are shown b y po ints, 

results by Hamer et al.l^are represented by dash line, dash-dotted and dotted lines are drawn using GSWA formula (42 I and 
Lindgard's formula ( |41[ |, respectively. Spectra eip, e2p and esp are also shown, which are found in the first, the second and 
the third orders in JjD and which are given by Eqs. ( |13[ l, ( |31[ ) and ( |36| l, respectively. The critical value of the anisotropy 
D^^ ~ 5.65^^^-^ is marked, below which the phase with the antiferromagnetic long-range magnetic order is stable. 



the following set of parameters is used for interpretation of the majority of experimental data: 

D = 8.9K, 

J, = 2.2K, (45) 
J,j, = 0.18K, 

where Jz is the exchange constant along the chains and Jxy is the exchange coupling constants between chains inside 
one tetragonal sublattice. Interaction between tetragonal sublattices is supposed to be negligibly smalP^ in the most 
of considerations. 
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FIG. 6. (Color online.) The ground state energy of the model |T| on the square lattice with 5 = 1 and antiferromagnetic 
exchange between only nearest neighbors. Monte-Carlo data of Refs.'^^'^ are shown by points, results by Hamer et alF^ are 
presented by dash line, the solid line is drawn using Eq. (37 1, which is obtained in the third order in J/D. 



DTN behavior attracts special attention near two QCPs in magnetic field H applied along the hard z axis. The 
first QCP H = Hci separates the paramagnetic and the canted antiferromagnetic phases and the second one separates 
the canted antiferromagnetic phase and the ferromagnetic one in which all spins are parallel to the field. Values of 
these critical fields are defined by the following exact relations (see, e.g., Refs.'^^'^Ii); 

Hcl = Epo , 



H,2^2SJo + Vo + Di2S~l), 



(46) 
(47) 



where Cp is the spectrum at = and V is the exchange coupling between sublattices considered below which has 
been neglected so far. The following values of the critical fields are obtained experimentally in Rei; ^^ ". ^^ ' 



2.05 T, 



ttDTN 
J^cl 

Hg^^ = 12.175 T. 



(48) 



In the present paper we focus on analysis of the DTN elementary excitation spectrum at zero magnetic field, which is 
observed in neutron experimeniP^ and shown in Fig. [7|^a)- (c). It s current theoretical interpretation looks inconsistent. 
Eqs. (|42])-(|44]) are used for the spectrum analysis in Refsl^SSH. As is shown in Ref.'^, Eqs. (|42])-(|44| describe DTN 
spectrum very well with parameters D = 8.12 K, J^y = 0.17 K and = 1.74 K which differ from those used in the 
literature now ( [45| (see the dash line in Fig. 7]ja)-(c)). However in the subsequent paper^*' these parameters were 
declined because the value of the critical field Hc2 found using Eq. ( 47 ) with these parameters differs significantly 
from the experimentally obtained value. The set of parameters (45) is proposed in Ref.^°, which has been used up to 
the present. In particular, in the recent papeil^S Eqs. (42|-(44) are used with the conventional parameters (45 1 for 
the spectrum analysis in the paramagnetic phase near the antiferromagnetic vector pp. However, as it is seen from 
Fig. [7]ja)-(c) (the dash-dotted line), GSWA with this set of parameters describes unsatisfactorily the spectrum of 
short-wavelength excitations. 

Because the body-centered tetragonal magnetic lattice of DTN is a Bravais lattice, all expressions obtained above 
are applicable to DTN. The spectrum obtained using Eq. (36) with the conventional parameters (45) for DTN is 
presented by the solid line in Fig. [7]^a)-(c). It is seen that the agreement with experimental data is poor. At the same 
time the method of the spectrum calculation proposed in the present paper looks suitable for DTN. It is illustrated by 
Fig.[7]jd)--(f), in which the second and the third order corrections in J/D to the spectrum and the first order spectrum 
eip given by Eq. ( 13 ) are presented. It is seen that the third order corrections are 3-5 times smaller than the second 
order ones in almost the whole Brillouin zone except for the vicinity of the antiferromagnetic vector, where they are 
almost equal but still remain much smaller than eip. 

Thus, the theoretical description of DTN needs revision that is indicated also by recent ESR-experimental data. As 
it is mentioned above, the interaction between spins from different DTN sublattices is usually ignored. But the data 
of the recent ESR experiment?^ indicate that it should be taken into consideration: one of the spectrum branches 
has a gap in the canted antiferromagnetic phase at Hd < H < Hc2 and the optical mode is slightly split. At the 
same time the model (Tj) without the inter-sublattice interaction has the doubly degenerate spectrum (due to two 
equivalent magnetic sublattices). 
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FIG. 7. (Color online.) (a)-(c) Elementary excitation spectrum of DTN along three directions in the Brillouin zone at zero 
field. The data of the neutron experimenlp^ at T = 80 mK are shown by points, dash and dash-dotted lines are drawn using 
Eqs. (|42|)-p4| with parameters D = 8.12 K, J^y = 0.17 K and Jz = 1.74 K (which are proposed in Ref.^^l) and with the 



conventional set of parameters (45 1, respectively, solid lines are drawn using Eq. (361 with the conventional parameters (45 1. 
(d)-(f) The spectrum of the first order in J/D given by Eq. ( 13 1 and the second and the third order corrections in J/D to the 
spectrum found using Eqs. (31 1 and (361 and the conventional set of parameters (451. 



In our previous paper'*' the spin- wave approach and the magnon Bose-condensation technique (near -ffc2) are used 
for analysis of the canted antiferromagnetic phase. It is shown that the inter-sublattice interaction of the form 



(49) 



where index i labels sites of one sublattice and j labels sites of another sublattice nearest to i, leads to the effects 
observed in the ESR experiment (the gap in one of the spectrum branches and the optical mode splitting). Unfortu- 
nately the lack of experimental data near Hc2, large anisotropy and quasi- ID nature of DTN did not allow to make 
reliable quantitative predictions about the value of Vij. It is just shown in Ref.'^that if Vij = V, then ^ ~ 0.1 K. 
In the present paper we continue our discussion started in Ref.^S and fit the elementary excitation spectrum using 



the least square method by varying parameters D, Jz, Jxy and V, while keeping the critical fields given by Eqs. (46) 



and (47 1 to be equal to the experimentally obtained values (48). As a result of this fit we obtain the following set of 



parameters that differs noticeably from the conventional one ( 45 ) 



D ^ 7.72 K, 
Jz = 1.86 K, 
Jxy = 0.2 K, 
V = 0.1 K. 



(50) 



The spectrum obtained is presented in Fig. [s] that is in good agreement with the neutron data of Ref.^l^. The inter- 
sublattice interaction removes the double degeneracy of the spectrum in DTN. This splitting is zero in (a) and (b) 
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panels of Fig.[8j but it is clearly seen in panel (c). It should be noted that the upper branch of the spectrum in Fig.|8][c) 
goes beyond the experimental error near p = 0. However, the branch splitting value Sep = e+ — Cp is very small 
compared to €p+€~ and it appears to be given in DTN by a slowly convergent series in J/D: Sep — (1.7— 1.0 + 0.4) K, 
where the first, the second and the third terms stand for the values of corresponding corrections in J/D. Then, one 
has to take into account higher order terms in J/D in order to find the small value of Sep in DTN that is out of the 
scope of the present paper. 

Note also that the value oi V ~ 0.1 K found above agrees with the estimation V ~ 0.1 K obtained as a result of 
our consideratiorP^ of the canted antiferromagnetic phase. 



V. CONCLUSION 



We propose the new representation ([s]), Q for an integer spin S via bosonic operators, which is useful in describing 
the paramagnetic phase and transitions to magnetically ordered phases in magnetic systems with large single-ion 
easy-plane anisotropy. Using this representation, the diagram technique and treating the exchange interaction as a 
perturbation we obtain Eq. ( 36 1 for the elementary excitation spectrum of the model ([I]) in the paramagnetic phase in 
the third order of the perturbation theory (that is referred to as an expansion in terms of J/ D for shot) . Expression 
(37) is also found for the ground state energy in the third order in J/D. Eq. ( [36| coincides with that obtained in 
ReffSSl in the special case of a spin chain with S — 1 and the exchange interaction between nearest neighbors only. 
We recover Eq. ( 36 1 also at 5* = 1 using simpler spin representations ([3| , ([38]) and ([3| , (|39]) . 

Comparison with numerical results obtained in Refs.l^SIifl for the square lattice, S ~ 1 and the antiferromagnetic 
exchange between nearest neighbors only shows that Eqs. (36) and (37) work better than results of other analytical 
methods proposed so far. In particular, Eqs. (36) and (37) work very w ell when D is not very close to the critical 

.hase becomes stableP^SHZl (gee Figs. liHel). At D > D^ Eq. (|36| poorly 



value Dc, below which the antiferromagnetic pi: 
describes only the spectrum of long-wavelength quasiparticles. 

It is shown that Eq. ([36]) is applicable for the spectrum analysis of the intensively studied compound NiCl2- 
4SC(NH2)2 (DTN), which is described by the model ([T])"^^^. We show that Eq. ^ with the conventional set of 
parameters for DTN^° (45) describes the experimentally obtained spectruml^^ unsatisfactorily (see Fig. [t]). The new 
set of parameters ( |50] ) is proposed for DTN which provides a good description of the experimental spectrum (see 
Fig. [sj and reproduces the experimentally obtained critical fields valueJ^ ( 48 ) . In contrast to the conventional model 
proposed for DTN before we take into account also the exchange interaction (49) between DTN magnetic sublattices, 
which becomes apparent in the recent ESR experiments^. 

In the forthcoming paper we continue analysis of the model ([T]) using representations ([3]), ([4|, (38) and (39) and 
consider its behavior in the vicinity of QCP H — Hd at T 7^ 0. 
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